Deflection of suspended graphene by a transverse electric field 
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Abstract 

We investigate the electromechanical response of doubly clamped graphene nanoribbons to a 
transverse gate voltage. An analytical model is developed to predict the field-induced deformation 
of graphene nanoribbons as a function of field intensity and graphene geometry. This model is 
validated thought atomistic simulations using the combination of a constitutive charge-dipole model 
and a pseudo-chemical many-body potential. As a newly observed effect of electric polarization, 
this field-induced deflection allows the graphene to oscillate at its natural frequency, which is found 
to decrease dramatically with increasing graphene size. 
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If we bring a glass rod electrically charged by rubbing with silk near a hair, the hair will 
be attracted to the rod. Here we demonstrate a similar electrostatic effect occurring in a 
graphene nanoribbon, a one-atom-thick carbon crystalline layer, which has been shown to 
have interesting electronic properties [1, 2j at nanoscale. In particular, graphene's electronic 
gap tunable in external electromagnetic fields makes it desirable for numerous applications 
in nanoelectronic devices [3llll[5l[6l[71[8l[9]. In such devices, the graphene is usually 
suspended between supports and exposed to an external electric field [10]. Since electronic 
transport properties of graphene are highly sensitive to the change of atomic structure [U |TT] , 
understanding its mechanical behaviors in response to an applied electric field is a crucial 
part for the development of graphene-based electronic and electromechanical devices. 

In recent experiments, vibrations induced by an electric field were observed in carbon 
nanotubes (CNTs) [T2]. This property has then been exploited in the design of a number 
of CNT-based nanoelectromechanical devices [I3j. It would not be surprising if graphene 
exhibits similar structural instability in an electric field as nanotubes do [H], in view of 
the large similarity in their in-plane electric polarizabilities [12] and mechanical properties 
[16] . Indeed, as an extremely thin membrane with well-defined electronic properties and 
low lateral stiffness, a suspended graphene sheet should be an ideal material for use in an 
electromechanical oscillator or resonator. 

When an electric field is applied across a thin neutral molecule, a moment of force acting 
on the molecule will result of electric polarization. It bends the molecule to the 

field direction [T7] . By such a deflection, electrical potential energy is converted and stock in 
the molecular structure. The mechanism of this energy conversion relies on the interaction 
between the electric field and the polarized charges on the molecule. This energy can be 
released in a form of mechanical oscillation by removing the field as far as the deflection of 
the molecule is elastic [18] . 

A natural question to ask at this point is, how does the shape of graphene change in 
response to applied transverse electric fields? In this work, we present an analytical approach 
to predict deformations of graphene nanoribbons (GNRs) by an electric field, demonstrating 
the coupling between the electric polarization and the mechanical resistance of graphene. 
This approach is validated through molecular simulations in view of the difficulty to establish 
an experimental quantification system of nanoscale electric polarization effects. Resonance 
frequency of GNRs with different length is measured by simulating vibration of GNRs by 
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means of molecular dynamics (MD). 

We performed molecular simulations [H] to compute the deflection of suspended GNRs 
by an electric field, minimizing the total potential energy of the system, which consists 
of two terms: an internal potential due to the C-C chemical bonds Uc, and an external 
potential Ue arising from the interaction with an external electric field. Uc is calculated 
using the adaptive interatomic reactive empirical bond order (AIREBO) potential function 
|19j . which has been used in recent studies on mechanical properties of CNTs [20] and GNRs 
[21] • Ue is computed using a constitutive Gaussian-regularized electrostatic model, in which 
each atom is modeled as an induced dipole p and a quantity of free charge q [22] , 

TV TV ^ N N 

i=l i=l i=l j=l 

N N ^ N N 

i=l j=l 1=1 j=l 

where N is the total number of atoms, x is the electron affinity, V and E stand for the 
external potential and electric field, respectively. T and T are the electrostatic interacting 
tensors. This model has recently been validated through electrostatic force microscopy 
experiments on CNTs [23] . Further details about the simulation can be found in Refs. [131 l2l] . 

Technically, no deformation will take place if an electric field is applied perfectly per- 
pendicular to the graphene surface, because the induced dipole is already parallel to the 
field direction hence the induced moment of electrostatic force is zero [25]. However, this 
ideal situation could be seldom attained due to previously predicted intrinsic height fiuctu- 
ations [26] and edge stresses [21] in graphene. During verification of these predictions using 
MD simulations based on the AIREBO potential, we have observed interesting elastic wave 
propagation on the surface of graphene at room temperature (see Fig. [l](b-d)). The typical 
speed and maximal vertical amplitude of this wave propagation in the graphene are found to 
be about 2 km/s and 0.2 — 0.4nm, respectively. The elastic waves (so called intrinsic ripples) 
make the shape of GNRs more or less naturally curved. This natural curvature provides 
possibility for an electric field to shift negative and positive charges to opposite directions 
in graphene. 

To show the physical principle of electric defiection, we depict in Fig. [T] (a) the profile of 
polarized charges in suspended graphene by a transverse electric field (computational details 
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can be found elsewhere [2^). We see that positive and negative charges in a GNR are shifted 
to its top (center) and bottom (side), respectively. Two pairs of opposite forces F^iec arise 
from the electrostatic interaction between the field and polarized charges, and form a bending 
moment M acting on each half graphene. Mirror symmetry of the system indicates that the 
GNR can equally be deflected downward, depending on the initial curvature of graphene at 
the moment when the gate voltage is applied. 

The key to modeling this electromechanical behavior is to understand the force equilib- 
rium in graphene. We consider the force balance in a half GNR as shown in Fig. [2] (a), 
which highlights the membrane-like character of graphene's mechanical properties. The for- 
mation of a moment of electric force can be understood by results plotted in Fig. |2] (b). 
This force profile is quite different from a commonly used assumption of a uniform electric 
force distribution in simplified calculations [281 EH]. In general, the amplitude of driving 
electric force should be proportional to either the polarizability of GNR {Feiec oc a), or the 
square of field strength {F^iec oc E"^) [T3]. Since the longitudinal polarizability of graphene 
is usually much larger than the transverse one [50] («// >> «±), M created in each half of 
the graphene by a electric force distribution (Fig. |2] (b)) can be calculated analytically as 

M = E'^a/fsinecose (2) 

where 9 stands for the angle of deflection (Fig. [T] (a)). In our earlier work [31], it has been 
shown that a// of a GNR is roughly proportional to its width w or the square of its length 
when the graphene sheet is not too small. Thus, for the system that we study here, Eq. 
pleads to 

M = E^Bw{^f sine cose (3) 

where B = 0.069nm is a constant related to the dielectric constant of graphene. 

We now estimate the moment of membrane force M* induced by the graphene shape 
change during the deflection. Considering the fact that graphene's stiffness in its atomic 
plane is about 30 times higher than that in the perpendicular direction [T6l 132] , the internal 
stress in a curved GNR can be supposed to be mainly due to its axial deformation e, which 
can be approximated as : e = a/Y ^ l/cos{e) — 1. Thus, M* can be expressed as 
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M* ^ -AY-tan9(^—-l) (4) 
2 ^cos^ ' ^ ' 

where Y is Young's modulus and A is the cross sectional area. Here we use the value of their 
product AY 2120i(7 eV/nm recently measured from nanoindentation experiments [T6] . 
This value is in good agreement with that from theoretical calculations [21] for GNRs with 
either armchair or zigzag edges. With the aid of moment balance in graphene M = —M* 
(shown in Fig. [2] (a)), and by combining Eqs. |3]and|4| we obtain the governing equation of 
electrostatic deflection of suspended graphene as follows: 



C 1 1 

where C = 2AY/Bw = GOGbSeV/nm^ is a constant. Since w has been eliminated from this 
equation, we can first conclude that the deflection of a GNR is independent of its width. 
We note that increase of cos 6 with decreasing w can be expected for narrow GNRs due to 
the effects of depolarization and edge states. However, this variation can be neglected for 
the size of GNRs usually reported in experiments. This equation also suggests that, for a 
given deformation angle 6, required fleld strength E decreases linearly with increasing the 
square root of graphene length ^/L. It is important to note that, since the parameters of our 
charge-dipole model were fltted for metallic sp^ system, Eq. [s] may not be useful for very 
narrow graphene sheets (with the minimum lateral dimension < 6 nm) which have large 
band gaps due to edge states ^ 

We now consider a particular case of small deflection (typically 6 < vr/12), for mainly 
covering experimentally reported resonance amplitude [101 EH]- In such a case Eq. [s] implies 
that 6 is roughly proportional to E since geometric simpliflcation gives cos~^ 6 — cos~^ 6 = 
tan (0.56*) tan cos^^ The maximal deflection 6 (Fig. [l](a)) can be approximately calcu- 
lated as: S = L tan 9/2. Thus, the case of small deflection from Eqis] leads to 



The analytical model presented in Eq. |6] is validated through series of molecular simula- 
tions for GNRs of different lengths. Results shown in FigjS] remark a quantitative agreement 
between analytical prediction and simulation data. It can be seen that 6 roughly follows 
a linear relationship with E and rapidly increases with L, as predicted by Eq. [Gj The 
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slight difference between tlie slopes of the curves is supposed to be mainly due to geometric 
approximations and the difference between experimentally measured Young's modulus and 
the one predicted from simulations. 

We now estimate the possibility of graphene's free oscillation induced by this electrostatic 
deflection by means of MD. In such simulations we apply a transverse electric field to deflect 
a suspended GNR, then the field is removed for inducing free oscillation of the GNR around 
its equilibrium position [36] . We simulated the vibration of graphene during a period of time 
(some ns), and we counted the number of oscillation from the output of the simulations. We 
found that the GNR oscillates with different harmonics, depending on its length. Only short 
GNRs {L < 12 nm) were observed to oscillate at fundamental harmonic and other GNRs 
oscillate at high-order ones. This is because the field-induced resonance in a GNR is a 
combination of transverse and longitudinal waves (see MD simulation video-recording [31]), 
while the longitudinal wave has no space to propagate in very short graphene hence only 
the fundamental harmonic was observed. We note that the vibration mode and frequency 
can also vary with the initial tension from supports due to different nature of fixation [T5] . 

We measured the natural frequency / of GNRs with different sizes. From Fig. |4]we can 
see that / decreases rapidly with increasing L. This length dependence is comparable to 
those found for single- walled CNTs [35]. In general, graphene's free oscillation is tunable by 
applying an alternating voltage, by which the mechanical vibration can be either enhanced 
or attenuated in terms of adjusting the AC frequency [M]- Theoretically speaking, the 
suspended structure could even be destroyed if the gate frequency is close to the natural 
frequency of the graphene. 

In conclusion, we have demonstrated that a suspended GNR can be deflected by applying 
a transverse electric field. The strong correlation between the field strength, graphene size 
and induced deflection has been concluded in an analytical model, which is validated via 
simulations. It was found that the deflection of GNRs is roughly proportional to the field 
strength for small deformations, and increases with the graphene length. The graphene's 
resonance frequency induced by this electrostatic deflection is found to decrease dramati- 
cally with the graphene length. These results suggest new potential applications of GNRs 
in electromechanical resonators or oscillators, which allow direct conversion from electric 
potential energy to mechanical energy in nanoscale. 
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FIGURES 




FIG. 1: (Color online) (a) Topographic diagram of charge distribution on a suspended GNR 
{L = 20nm) in a transverse electric field E from molecular simulations. The color scale of atoms 
is proportional to the induced-charge density. .Fekc ^-nd — F'eiec stand for the electrostatic force 
arising from the interaction between the charges and the field, which make the GNR deflected. 
6 and 6 stand for the maximal amplitude and angle of deflection, respectively. The electric field 
can be generated by applying a gate voltage between parallel capacitors, (b-d) Representative 
atomic configurations of the graphene at room temperature from MD simulations before the field 
is applied. The dashed lines represent the cross section shapes and the arrows show the direction 
of elastic wave propagation. 
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FIG. 2: (Color online) (a) Schematic of force/moment balance on a suspended GNR {L/2 ~ lOnm, 
from the left end to middle) in a transverse electric field E. a A stands for a force pair of internal 
stress arising from the deflection of graphene. M and M* correspond respectively to the moment 
of electric force and internal stress, (b) Profile of electric force acting on atoms in different positions 
along a half of the graphene with different field strengths E. 
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FIG. 3: 5 vs E for GNRs with different lengths L. The symbols represent simulation data 
and the lines stand for those predicted by analytical model using Eqj5] (unit conversion: IV/nm 



0.833 v^ey/nm3). 
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FIG. 4: Resonance frequency / of GNRs of different lengths L. Inset shows that GNRs oscillate 
with different harmonics modes. The circles represent simulation data and the curve stands for a 
best-fitting equation: / = -6.532 + 778. IL'^ - 482.4L-2. 
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